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SUMMARY 


V  The  statistical  variation  of  load  and  strength  is  described  by  a  three 
parameter  Weibull  distribution.  The  Weibull  parameters  are  evaluated  by  a  least 
square  analysis  and  a  method  is  presented  which  allows  confidence  bounds  to  be 
assigned  to  these  quantities.  A  Monte  Carlo  analysis  is  used  to  calculate  the 
reliability  of  the  structure  from  the  load  and  strength  distributions. 
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INTRODUCTION 


The  Introduction  of  new  materials  for  load  carrying  applications  often 
presents  the  designer  with  problems  in  assessing  the  reliability  of  the 
structure.  Although  considerable  advances  have  been  made  by  material  scientists 
in  the  development  of  new  materials  capable  of  operating  in  severe  environments, 
it  is  only  during  recent  years  that  engineers  have  considered  introducing  probe-' 
bilistic  methods  of  design  into  structural  analysis. 

With  the  advent  of  engineering  plastics,  fibre  reinforced  composites, 
graphites,  ceramics  and  high  strength  steels,  many  materials  are  being  developed 
which  have  inherent  variability  in  thdir  properties.  For  example,  most  graphites 
and  ceramics  have  random  distributions  of  flaws,  both  through  their  volume  and 
on  their  surface,  which  lead  to  significant  strength  variations  between  nominally 
identical  components  subjected  to  identical  states  of  stress^.  It  is  a  rare 
event  in  the  aerospace  industry  to  find  a  material  property  or  an  applied  load 
which  is  not  subject  to  some  significant  statistical  variation. 

The  relative  strengths  of  most  structures  vary  statistically  and  direct 
measurement  of  the  strength  of  each  member  could  be  uneconomic  and  impracticable. 
To  overcome  these  difficulties  the  structural  reliability  may  be  assessed  from  a 
mathematical  analysis  based  on  the  test  results  from  a  small  sample  of  components. 

Once  the  statistical  variations  of  strength  and  load  have  been  approximately 
modelled  by  some  distribution  function,  the  re 1 i at^i 1 1 ty  of  the  structure  may  be 
evaluated  by  analysing  the  interaction  between  the  two  distributions.  At  the 
expense  of  imposing  certain  penalties  reliability  may  be  improved  by  reducing 
the  applied  load,  by  raising  the  mean  strength  of  the  material  substantially,  or 
by  reducing  the  variability  of  the  material  under  development.  If  the  material 
is  made  for  a  particular  Service  application  strength  can  only  be  increased  by 
thickening  cross  sections.  This  is  not  only  expensive  but  may  increase  the 
weight  to  an  unacceptable  level  in  certain  aerospace  applications.  Material 
variability  can  be  reduced  if  there  is  a  choice  of  materials  available.  The 
materia]  variability  of  some  composites  can  be  reduced  without  adversely 
affecting  the  strength  and  other  properties,  by  changing  the  material  compo¬ 
sition  slightly. 

The  accuracy  of  predicting  the  extreme  values  of  a  statistical  variable 
strength  and  load  depends  on  the  fidelity  of  the  mode),  and  on  the  amount  of 
data  available  and  their  integrity.  Predictions  based  on  the  analysis  of 
several  hundred  tests  may  be  expected  to  give  reliable  results  but  those 
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evaluated  from  a  few  tests  are  Inevitably  less  precise.  In  these  circumstances 
predictions  are  made  using  conservative  estimates  for  the  population  parameters. 

Hany  statistical  distributions  are  commonly  used  to  describe  experimental 
data  and  the  choice  of  distribution  in  the  field  of  reliability  is  by  no  means 
unique.  The  normal  distribution  might  appear  to  be  suitable  and  has  the  advan¬ 
tage  that  the  evaluation  of  the  confidence  limits  for  small  sample  data  is  well 
documented2.  However,  the  normal  distribution  is  symmetric  and  is  inappropriate 
for  fracture  data  with  skewed  distributions.  It  is  also  double  sided,  extending 
to  infinity  in  both  positive  and  negative  directions,  which  implies  a  finite 
probability  value  when  the  distribution  variate  is  zero.  In  a  fracture  analysis 
application  the  normal  distribution  would  yield  a  finite  probability  value  for 
zero  load  which  is  clearly  physically  inadmissible,  but  it  is  negligible  in 
practical  terms  and  allows  the  normal  distribution  to  be  used  in  many  cases. 

A  distribution  which  is  bounded  to  the  left  and  allows  for  positive  and 
negative  skewness  Is  the  one  due  to  Weibull  ,  which  has  wide  application  in  the 
analysis  of  fracture  and  material  strength  and  is  adopted  in  this  analysis.  A 
Weibull  distribution  is  also  chosen  to  describe  the  statistical  variation  of 
loads.  Minor  modifications  would  be  needed  if  the  load  variations  were  described 
by  some  other  distribution  function. 

In  Section  2  of  this  report  the  strength  and  load  variations  are  each 
described  by  a  three  parameter  Weibull  distribution  function.  A  least  square 
analysis  is  presented  which  allows  the  various  Weibull  parameters  to  be  calcu¬ 
lated  from  the  known  sample  data.  Confidence  limits  are  assigned  to  the  para¬ 
meters  through  an  analysis  outlined  in  Section  3*  In  Section  A  a  Monte  Carlo 
method  is  presented  which  allows  the  reliability  of  a  structure  to  be  derived 
from  the  strength  and  load  distributions. 

2  WEIBULL  ANALYSIS 

When  the  strength  capability  of  a  structure  is  exceeded  by  an  applied  load 
the  structure  fails  and  no  longer  meets  its  original  design  requirement.  If  the 
variability  associated  with  strength  and  load  is  small  it  is  possible,  to  a 
reasonable  degree  of  accuracy,  to  express  failure  by  an  Inequality  Y>X.  Here  X 
and  Y  denote  strength  and  applied  load  respectively.  There  are  many  situations 
where  X  and  Y  are  not  deterministic  quantities  but  vary  statistically  accord¬ 
ing  to  some  distribution.  In  this  case  a  criterion  of  failure  expressed  in  terms 
of  some  inequality  becomes  inadmissible  and  must  be  replaced  by  an  expression  for 
probability  of  failure. 


The  Weibull  distribution^  is  extensively  used  in  engineering  and  design, 
and  its  application  to  the  strength  analysis  of  brittle  materials,  composites 
and  polymers  is  well  documented1*’^'^'^'^.  If  the  three  parameter  equation  is 
used  an  allowance  may  be  made  for  threshold  value,  skewness  and  scale.  The 
threshold  value  defines  a  point  In  the  range  below  which  the  probability  of  a 
critical  event  is  zero.  Above  the  threshold  value  the  probability  of  a  critical 
event  conforms  to  the  Weibull  distribution. 

in  this  analysis  it  is  assumed  that  both  the  strength  and  applied  loads  are 
distributed  according  to  the  three  parameter  Weibull  distribution.  It  is  antici¬ 
pated  that  in  most  applications  this  assumption  will  be  reasonable  for  the 
strength  distribution  but  there  may  be  certain  circumstances  in  which  a  different 
distribution  is  required  to  model  the  statistical  variation  of  the  load.  In  this 
event  the  analysis  will  require  minor  modification. 

The  cumulative  probability  functions  for  the  applied  load  and  strength  are 
given  by  the  respective  relationships 


In  the  above  equations  P^  and  Py  are  the  cumulative  probability  values,  M^ 
and  My  are  the  Weibull  moduli,  XQ  and  Yg  are  the  normalising  factors,  Xy 
and  Yy  are  the  threshold  values  of  the  respective  variates  X  and  Y  .  The 
Weibull  moduli  and  normalising  factors  respectively  govern  the  skewness  and 
scale  of  the  distribution.  The  threshold  values  are  responsible  for  origin 
shift.  In  terms  of  these  Weibull  parameters  the  distribution  means  and  coeffici¬ 
ents  of  variation  are  given  respectively  by^ 
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X 

Y 


XT  +  X0  r  0  +  1^  ] 

yt  +  Yo  r  (1  +  7%  ) 


(3) 


The  Veibull  parameters  appearing  in  equations  (1)  and  (2)  are  population 
values.  For  most  engineering  applications  these  quantities  are  not  available 

but  will  need  to  be  estimated  from  small  sample  data  to  a  given  confidence  level. 
If  the  small  sample  data 


X  -  Xi  *  Xi+1  >  X, 


Y  *  Y. 


*  vi+i  51  Y1 


V  i**1 ,2 . N  , 


(5) 


are  available,  where  for  convenience  the  size  of  each  sample  Is  taken  to  be  N  , 
the  Weibull  parameters  may  be  evaluated  by  a  least  square  analysis10.  This  is  ’ 
achieved  by  minimising  the  error  terms 


6X(,)  -  X  (Px(i))  -  X.  -  X,  -  X, 
«y(i>  -  V  (PV0>,  -  Y,  -  i,  -  », 


> 


(6) 


with  respect  to  the  unknown  Weibull  parameters.  In  equation  (6)  X.  and 
Yj  satisfy  the  transformed  Weibull  equations 
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1 


Computations  have  shown  that  the  assigned  experimental  probability  values 
provided  by  equation  (8)  lead  to  unbiased  estimates  of  the  Wei  bull  parameters^. 

It  can  be  shown  that  the  minimisation  of  equation  (6)  Is  equivalent  to 
solving  the  least  square  equations 
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Equations  (9)  and  ( 10)  may  be  solved  simultaneously  to  express  Xq  and 
XT  in  terms  of  the  Weibull  modulus  Mx  .  Omitting  the  details 


Substitution  for  Xq  and  X-j.  into  equation  (11)  leads  to  a  transcendental 
equation  for  .  Although  its  final  form  is  too  complex  to  be  written  in 

this  text  the  resulting  equation  may  be  readily  solved  numerically  on  a  high 
speed  computer.  In  the  root  solving  routine  adopted  a  starting  value  of 
»  1.0001  was  chosen  because  for  M  *  1  the  Weibull  equation  degenerated 
into  a  physically  unrealistic  form  for  the  studies  under  investigation.  The 
solution  for  the  other  parameter  set  My  ,  Yq  ,  Yy  may  be  obtained 
simi larly. 
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3  CONFIDENCE  LIMITS 


In  equations  (1)  and  (2)  the  Weibull  parameters  refer  to  population  values. 
The  Weibull  parameters  calculated  from  equations  (9)  to  (15)  are  small  sample 
estimates  and  only  approach  the  population  values  as  the  size  of  the  sample 
increases  to  infinity.  In  practice  the  sample  size  may  be  small  and  the  Weibull 
parameters  evaluated  from  this  limited  data  may  differ  significantly  from  the 
population  values,  in  these  circumstances  it  is  useful  to  be  able  to  define, 
for  a  given  level  of  confidence,  upper  and  lower  limits  of  the  deviation  of 
these  parameters  from  the  population  values. 

When  all  three  Weibull  parameters  vary  statistically  the  estimation  of  the 
confidence  bounds  requires  individual  attention  for  each  new  set  of  parameters 
encountered.  If  the  population  threshold  values  are  known  or  may  be  prescribed 
considerable  simplifications  occur.  The  confidence  limits  for  the  remaining 
two  Weibull  parameters  may  be  evaluated  from  a  Monte  Carlo  analysis  of  a  single 
system  of  equations  and  these  results  hold  universally  for  all  Weibull  distri¬ 
butions. 

Before  proceeding  further  it  is  important  to  note  that  the  least  square 
equations  may  be  developed  from  any  function  of  the  differences  defined  by 
equation  (6).  In  particular  the  quantities 


(18) 


may  be  used  to  carry  out  the  least  square  analysis.  With  these  assumptions 
explicit  expressions  for  Mx  »  X0  »  *Y  *  *0  may  be  derived  ,n  terns  of 
the  threshold  values  and  Y^  .  The  threshold  values  are  found  as  roots 

of  two  transcendental  equations.  Unfortunately  these  values  are  highly  sensitive 
to  small  changes  in  data  and  in  many  cases  the  calculated  Weibull  parameters  are 
spurious.  If  the  quantities  and  Y^  are  prescribed,  an  analysis  based  on 
equation  (18)  yields  reliable  and  accurate  results  for  all  cases. 

If  the  threshold  value  Xy  is  known  a  least  square  analysis  based  on 
equation  (18)  yields  the  parameters 


M. 


(S) 


(19) 


*o(S’  ■  ««P  i 


N  N  w  n 

Vj2yE  ,n  (W  ‘  2  V|  V|  ,n  (xi~xt 
“mt  frr  .... 


N  N 

I -I 

1=1  1=1 

N  N 


i=1  1=1  1=1 


(20) 


where 


In  Jn^  1/(1-Px(i))  j 


(21) 


In  the  above  equations  the  small  sample  values  »  X(/^  are  introc,uced 
to  distinguish  them  from  the  population  values  ,  Xp  Making  the  change 
of  variable 


M 


X 


reduces  the  Weibult  distribution  to  the  exponential  distribution 

Px  =  1  -  exp  (-x)  . 


Substi tut ing  equation  (22)  into  equation  (19)  leads  to 


(22) 


(23) 
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(24) 


N  N  N 


The  quantities  x.  are  derived  from  equation  (22)  with  X  set  equal  to 

X.  i <*1 ,2 . N.  If  the  terms  in  the  brackets  appearing  in  equation  (24)  are 

expanded, 


N  N  N 


When  equation  (22)  is  substituted  into  equation  (20)  and  the  resulting  expression 
simplified  it  can  be  shown  after  some  manipulation  that 


In  the  above  two  equations  the  right  hand  sides  in  each  case  are  independent 
of  the  Weibull  parameters.  They  depend  only  on  the  distribution  of  x  ,  the 
random  variable  of  the  exponential  distribution.  This  implies  that  a  Monte  Carlo 
analysis  applied  to  equations  (25)  and  (26),  for  various  values  of  N  ,  will 

M  (S) 

result  in  distributions  of  ancJ  (X0^/Xq)  *  which  are  universally 

valid  for  all  Weibull  parameters. 

Random  samples  of  any  size  N  may  be  generated  by  evaluating  the  inverse 
of  the  exponential  distribution 
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X.  -  -  In  (1-P  ' 

I  x 


1-1,2 . N 


for  a  set  of  random  numbers  P 


between  zero  and  unity.  These  quantities 


may  then  be  substituted  into  equations (25)  and  (26)  and  the  values  of 

M  (S) 


Mx(S)/Mx  and  X 


calculated  for  that  particular  sample.  The 


quantities  V.  are  known  and  given  by  equations  (8)  and  (21).  If  this  proce¬ 
dure  is  repeated  many  times  for  the  same  sample  size  N  ,  a  typical  number  of 

m  (S> 

fsl  /el  ny 

times  being  2500,  distributions  of  M^'  '/M^  and  (XqV  /Xq)  may  be 

established.  It  is  then  a  trivial  matter  to  evaluate  the  confidence  limits 
directly  from  the  percentiles  of  the  ranked  distributions.  A  Monte  Carlo 
diagram  showing  this  sequence  is  given  in  Fig.  1. 

The  above  procedure  may  be  repeated  for  any  sample  size  N  .  Consider¬ 
able  savings  in  computer  time  may  be  achieved  if  the  results  can  be  expressed 
as  some  equation  of  N  and  the  confidence  bounds.  To  establish  such  relation- 

M  <S> 

(c\  (c)  “v 

ships  500  values  of  M^'  and  (Xq  /Xq)  were  obtained  10  times  for 

each  value  of  N  ranging  from  10  to  100  in  steps  of  10.  The  first  three 
moments  of  the  distribution  of  results  were  evaluated  to  measure  the  mean, 
variance  and  skew.  An  examination  of  these  parameters  and  the  ordered  results 
indicate  a  log-normal  distribution  in  both  cases.  The  variables  were  converted 
to  their  logarithms  and  subsequent  analysis  confirmed  the  distribution  to  be 
log-normal.  The  average  values  of  the  mean,  variances  and  skew  moments  are 
listed  in  Tables  1  and  2. 

In  Table  1  it  is  shown  that  for  each  sample  size  both  the  mean  and  variance 

are  approximately  equal  to  the  reciprocal  of  the  sample  size.  The  skew  moments 

are  small  and  appear  to  be  distributed  about  zero  indicating  symmetrical  distri- 

(S) 

buttons.  The  distribution  of  !n(M^'  VM^)  can  be  approximated  very  closely 
by  a  normal  distribution  with  mean  and  variance  equal  to  l/N  implying  that 
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In 


(28) 


y  (S) 

JL _  .  1  *  7  1 

l  m  J  I  4T 


where  Z  is  the  normal  standard  deviate  corresponding  to  the  prescribed  value 
of  the  confidence  limit.  Confirmation  of  this  result  was  obtained  by  a  Konte 
Carlo  analysis  in  which  the  values  of  My^/H)(  were  evaluated  for  various 
values  of  N  from  a  particular  Weibull  distribution  whose  population  parameters 
were  ■  2  ,  Xg  ■  1  and  Xy  •  0  .  The  agreement  between  the  two  methods  is 

very  close,  as  shown  in  Table  3. 

m  <S) 

(S)  A 

The  mean  value  of  the  ratio  In  (Xg'  /Xg)  does  not  appear  to  be 

related  to  the  sample  size,  as  Table  2  shows,  but  being  smaller  than  its  standard 
deviation  it  may  be  assumed  to  be  zero.  The  standard  deviation  of  the  ratio  has 
a  log-linear  relationship  with  the  sample  size  N  and  may  be  closely  approxi¬ 
mated  by  the  equation 


with  a  correlation  coefficient  of  0.9987.  In  the  above  equation  oy  denotes 

M  <S>  ° 

(S)  A 

the  standard  deviation  of  ln(Xg  /Xg)  .  It  therefore  follows  that 


ln(X0<S)/X0)  X 


-  |  N-  ™  Z 


where  Z  Is  the  normal  standard  deviate  corresponding  to  the  prescribed 
confidence  limit  value.  A  Monte  Carlo  analysis  gave  similar  results,  as  shown 
in  Table  4. 

4  RELIABILITY  EVALUATION 

Equation  (7)  allows  the  distribution  variates  X  and  Y  to  be  expressed 
in  terms  of  the  Weibull  population  parameters.  In  practice  these  quantities  are 
not  usually  known  and  instead  only  estimates  from  small  sample  values  are  avail¬ 
able.  Provided  that  the  threshold  values  XT  and  Yj  can  be  prescribed,  the 
population  values  for  the  remaining  Weibull  parameters  lie  within  the  confidence 
intervals 
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Mv(s)  exp 


Xq(S)  exp 


3  H_Vf7?  7 
^  ' 


c  X0  <  X0(S)  exp 


-  vT73 


with  a  similar  pair  of  equations  for  My  and  Yg  .  The  quantities  and 

Z ^  are  the  normal  standard  deviates  for  the  lower  and  upper  confidence  limits 
respectively. 

Failure  occurs  when  the  applied  loads  exceed  the  component  strength,  and 
the  structural  reliability  may  be  obtained  by  considering  the  distribution  of 
their  differences.  A  conservative  estimate  of  failure  probability  may  be 
obtained  by  taking  upper  and  lower  confidence  estimates  of  the  population  values 
of  the  respective  load  and  strength  distributions.  The  minimum  and  maximum 
values  of  the  respective  variates  X  and  Y  follow  from  substitution  of  the 
values  of  the  inequalities  given  by  equations  (30  and  (32)  into  equation  (7). 
Omitting  the  details  the  minimum  and  maximum  values  of  X  and  Y  are  given  by 


(min)  a  ,/MX 
XT  +  X0  X 


(min) 


Y  -  yt  +  yo 


(max)  „  ,'»V 


(min) 


where 


(min)  y  (S) 

''A 


(<~«>  .  „o<S>  txp  |  _3  ™  2jJ  , 


\ 

\ 


(37) 


«x  •  '"(  ,/(,-px>  )  • 

(»Y  -  ln|  1/(1-PY)  |  , 


(38) 

(39) 

(40) 


Since  the  population  parameters  are  constant  only  one  confidence  limit  value 
may  be  chosen  for  each  parameter.  The  values  chosen  are  those  which  minimise 
the  strength  variable  In  the  lower  probability  region  (i.e.  <  1)  and  maximise 

the  load  variable  in  the  higher  probability  region  (i.e.  <*Y  >  1).  The  choice  of 
values  for  and  aY  Is  reflected  in  the  derivation  of  the  above  equations. 

The  probability  that  Y  exceeds  X  can  now  be  evaluated  by  a  Monte  Carlo 

analysis  in  which  the  values  of  Y  and  X  are  selected  by  generating  random 

values  for  and  PY  The  number  of  occasions  that  Y  exceeds  X  , 

expressed  as  a  ratio  of  the  total  number  of  trials,  gives  the  required  failure 

probability.  If  the  anticipated  probability  values  are  high  the  method  described 

above  may  be  acceptable,  but  for  low  probability  values  the  method  is  both 

inefficient  and  uneconomical  and  an  anticipated  probability  value  of  10  ^  would 

_g 

require  10  trials  to  give  moderate  results. 

An  alternative  method  is  to  analyse  the  distribution  of  the  differences 
defined  by 


0,  -  Y}  -  X,  ,  1-1,2,  - N  .  (41) 

Confutations  have  shown  that  up  to  N  -  1000  the  assumption  that  the  distribution 
of  the  differences,  D,  ,  is  normal  is  not  rejected  by  the  Lilliefors  or 
Kolmogorov  Smirnov  tests  at  the  20%  level  of  significance  .  Under  these 
conditions  the  variate  0  can  be  expressed  In  the  form 

0  -  UD  *  2  oD  (42) 
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where  Wq  and  Op  respectively  denote  the  mean  and  standard  deviation  of  the 
difference  distribution  Dj  ,  1-1,2,  ....  N  .  The  quantity  Z  is  the  normal 
standard  deviate.  The  probability  that  Y  exceeds  X  is  found  from  equation 
(42)  by  first  setting  D  equal  to  zero  and  then  finding  the  probability  value 
corresponding  to  the  normal  deviate 


Z  -  -  -51 


(43) 


The  required  failure  probability  value  follows  from  the  relationship 

p['v>x:i  -  p£  0*0*1  -  1  -  P £ Z  -  -  •  (44) 

A  flow  diagram  describing  the  Monte  Carlo  simulation  for  calculating  P^^xj 
is  given  in  Fig.  2. 

To  assess  the  accuracy  of  the  above  procedures  population  parameters  were 
prescribed  on  the  assumption  that  the  two  distributions  for  X  and  Y  were 
normal.  The  population  parameters  were  chosen  to  give  a  high  failure  probability 
for  P^Y>x3  and  insignificant  probability  values  for  Pg  and  Py  as  X  and 
Y  approached  zero.  A  Monte  Carlo  analysis  was  used  to  generate  two  samples  of 
size  M  -  1000  from  the  population  parameters  characterising  the  above  normal 
distributions.  Each  sample  was  subsequently  divided  into  smaller  samples  of 
size  N  -  10,  20,  20,  50.  100,  200,  500.  The  failure  probability,  p[y>x3  , 
was  first  calculated  for  each  sample  size  N  to  a  specified  confidence  level  of 
952  by  the  non-central  1 11  statistic1^.  This  analysis  took  into  account  the 
total  interaction  between  the  two  distributions  Pg  and  Py  but  was  valid  only 
for  normal  distributions.  Wefbull  distributions  were  then  fitted  to  the  same 
sets  of  data  and  the  smell  sample  estimates  of  the  Weibull  parameters  calculated. 
Using  these  parameters  and  a  specified  Confidence  level  of  952  the  probability 
of  Y  exceeding  X  ,  P|~Y>xJ  ,  was  evaluated  by  the  Monte  Carlo  analysis  out¬ 
lined  at  the  beginning  of  this  section.  A  comparison  between  the  Monte  Carlo 
analysis  and  the  results  obtained  from  the  non-central  1 1*  statistic  are  shown 
in  Table  Agreement  is  close  in  all  cases,  the  Monte  Carlo  results  being 
slightly  more  conservative. 
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TABLE  1 


Distribution  parameters  of 


derived  from  a  Monte  Carlo  analysis 


Sample 

size 

Mean 

Variance 

Skew  moment 

10 

B 

BB 

BO 
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B 
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BB 
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BO 
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TABLE  2 


derived  from  a  Monte  Carlo  analysis 


Sample 

size 
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Variance 

Skew  moment 
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E 
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TABLE 


Distribution  parameters  of  I 


,  <j  (S) 

•(VI 


obtained  by 


Honte  Carlo  methods  from  a  Weibuli  distribution 
with  M  -  2  ,  XQ  -  1  ,  Xy  -  0  . 


Sample 


Hean  Variance 


4th  moment 


TABLE  4 


Distribution  parameters  of  In 


(S)  MX 


A  ~ 

j  obtained  by 


Honte  Carlo  methods  from  a  Weibuil  distribution 
with  H  •  2  ,  X.  •  I  .  i.  »  i  . 


Sample 

size 
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.0022 
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For  each  value  of  N  the  distributional  parameters  of  $00  samples  of  the 

u  ($) 

n  *ols!  ) *  i 

ratio  In  ■  J  were  obtained  10  times  and  averaged. 

Note,  for  a  normal  distribution  4.t-h.-lffg*£t  .  3.0  . 


variance 
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TABLE  5 


Probability  of  failure  calculations  for  various 
sample  sizes  at  a  confidence  level  of  95% 


Sample 

size 

Probability  of  failure  P Y> X ~ j 

Non  central  ' t * 
statistic 

Monte  Carlo 
analysis 

10 
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20 
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0. 1611 

0.1808 

200 

0. 1470 

mm 

500 

0.1358 

0.1347 

gg| 

0. 1276 

0. 1399 
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1  Monte  Carlo  diagram  for  confidence  limits 
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using  the  Monte  Carlo  similation 
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